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Localization properties of the doped Z2-topological insulator are studied by weak localization 
theory. The disordered Kane-Mele model for graphene is taken as a prototype, and analyzed with 
attention to effects of the topological mass term, inter- valley scattering, and the Rashba spin-orbit 
interaction. The known tendency of graphene to anti-localize in the absence of inter- valley scattering 
between K and K' points is naturally placed as the massless limit of Kane-Mele model. The latter is 
shown to have a unitary behavior even in the absence of magnetic field due to the topological mass 
term. When inter- valley scattering is introduced, the topological mass term leaves the system in the 
unitary class, whereas the ordinary mass term, which appears if A and B sublattices are inequivalent, 
turns the system to weak localization. The Rashba spin-orbit interaction in the presence of K-K' 
scattering drive the system to weak anti-localization in sharp contrast to the ideal graphene case. 



I. INTRODUCTION 

The concept of Z2 topological insulator was first intro- 
duced in a model for graphene [1], in the presence of both 
intrinsic and extrinsic (Rashba) spin-orbit interactions 
(called hereafter, Kane-Mele model) [H, The origin of 
Z2 symmetry lies naturally in the time reversal invariance 
of the underlying spin-orbit interactions, i.e., in the ex- 
istence of Kramers pairs. In the continuum limit, the in- 
trinsic spin-orbit interaction is represented by a so-called 
topological mass term (of size A, opening a spin-orbit gap 
2A), encoding quantized spin Hall effect. The latter oc- 
curs when Fermi energy is in the gap, and implies the ex- 
istence of a pair of counter-propagating gapless states at 
the sample boundary with opposite spins, often dubbed 
as helical edge modes [4]. The idea of "Z2" topological 
insulator stems from the observation that these helical 
edge modes are robust against weak perturbations, such 
as the extrinsic Rashba spin-orbit interaction (coupling 
strength: \r). Thus, non-trivial topological nature of a 
Z2 topological insulator is often attributed to the exis- 
tence of such edge modes, protected by Kramers degen- 
eracy. This paper, on the contrary, highlights its hulk 
property. Since real materials always have disorder, we 
investigate its transport property under doping using the 
framework of standard weak localization theory. 

Of course, the magnitude of spin-orbit interactions has 
always been questioned in graphene [i, 0, Q , leading to 
search for Z2 nature in a system of larger spin-orbit cou- 
pling y, . The existence of helical edge modes was first 
experimentally shown in a two-dimensional HgTe/CdTe 
heterostructure Q. Recall that in graphene two doubly 
degenerate Dirac cones appear at K- and K'- points in 
the first Brillouin zone [1], in contrast to a single pair 
of Dirac cones appearing at the F-point in HgTe/CdTe 
qauntum well. The first estimate of A and \r in the orig- 
inal paper of Kane and Mele: 2 A ~2.4 K, and Xr/2 ^ 
0.5 mK for a typical strength of perpendicular electric 
field E = 50 V/300 nm, provides a favorable condition 
for Z2 non-trivial phase [2]. This estimate was later 
shown to be too optimistic (for the occurrence of Z2 
phase) due to the specific geometry of s and p orbitals 



in graphene. According to Refs. [5], the estimated value 
of A {Xr) is much smaller (larger) than the original es- 
timation of Ref. [2]: 2 A -0.01 K, and Xr/2 - 0.13 K 
for the same electric field of £^ = 50 V/300 nm. On the 
other hand, a recent first-principle calculation suggests 
that d-orbitals play a dominant role in the gap opening 
at K and K' points [9]. As a result, the actual value of 
A might be somewhat intermediate between the previous 
estimates of [2] and [5, 6, 7], namely 2 A -0.28 K, Xr/2 - 
0.23 K per V/nm. The concept of Z2-topological insu- 
lator has also been extended to three space dimensions 
EO, [O, E E, [ii, [H, [H, [Til- a recent spin-ARPES 
study on Bi2Te3 reports on the experimental observation 
of a spin-helical two-dimensional surface state in such 
three-dimensional Z2-topological insulator fv^ . 

Focalization properties of the doped Kane-Mele Z2 in- 
sulator have been studied numerically [I81, lIB] • Ref. [iHl 
deduces a phase diagram in the {E, VF)-plane {E: energy, 
W: strength of disorder), in which a metallic domain ap- 
pears in valence and conduction bands with a finite width 
in As disorder is increased, these two extended do- 
mains in both bands approach to each other, and even- 
tually merge and disappear. A more subtle issue is the 
nature of the metallic state next to the Z2 insulating 
phase. It has been claimed flS*] that the system's Z2 sym- 
metry leads to an unconventional symmetry class. How- 
ever, an extensive study on the critical exponents 
has suggested that the weak anti-localization behavior of 
the doped Z2 insulator belongs to the conventional sym- 
plectic symmetry class. This paper addresses the basic 
mechanism how doped Z2 insulators acquire such unique 
localization properties. As a simple implementation of 
Z2-topological insulator, we consider Kane-Mele model, 
and in contrast to the numerical works of Refs. [isl, fl9|, 
we restrict our study to the limit of weak disorder. On 
the other hand, we pay much attention to the existence 
of valleys in graphene, since localization properties are 
much influenced by the presence or absence of scattering 
across different valleys in the Brillouin zone. The later is 
determined by the range of the impurity potential [g^. 

This paper is organized as follows. The Kane-Mele 
model is introduced in Sec. II. Then, we apply the stan- 
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dard diagrammatic approach to weak locahzation to the 
doped Kane-Mele model. In Sec. Ill we consider the 
case of vanishing Rashba SOI. Particular attention will 
be paid to different types of the mass term, (a) and (b), 
together with the presence/absence of K-K' scattering. 
Here, we will focus on unitary behaviors, which appear as 
a consequence of a finite lifetime acquired by Cooperons. 
Breaking or preserved effective time reversal symmetry 
will be the main issue of this section. Sec. IV is devoted 
to study on the effects of Rashba spin-orbit interaction. 
In the final section, we will summarize our results, and 
give interpretation to them in terms of the number of 
active species of effective spins [21i] . 



II. KANE-MELE MODEL 

The Kane-Mele model is given a status as a proto- 
type for various Z2 topological insulator models. It was 
introduced as a model for graphene in the presence of 
spin-orbit interactions [2|. The model is first defined on 
the hexagonal lattice, in the framework of tight-binding 
approximation. The continuum limit is then taken, in 
which the effective Hamiltonian becomes, 

Hkm = Hi -\- H^ + Hr 

Hi = hVpiPxCTxTz ^PyCTy) 

Ha = -Aa^TzSz 

Hr = -^{(JxTzSy - (JySj,), (1) 

where three Pauli's matrices, a, r , and s operate in dif- 
ferent spaces. Namely, <j acts on pseudo-spin specifying 
the A-B sublattices, r on the K-K' "valley spin" , and s 
on the real spin. Throughout this paper, we assume that 
A > 



A. Three story structure 

The Kane-Mele model defined as Eq. ([1]) has the fol- 
lowing three story structure (see TABLE I), correspond- 
ing to each term of Eq. ([1]): (i) graphene on its base, (ii) 
a topological mass term, encoding quantized spin Hall 
(QSH) effect, and finally (iii) the Rashba spin-orbit in- 
teraction \r. Let us first look into the role of these three 
fioors one by one. 



nearest neighbor hopping term in the tight binding ap- 
proximation. Such ideal (massless)_graphene shows weak 
anti-locazation behavior (2^, [2^, when inter- valley 
[K-K') scattering is irrelevant [23]. The absence or pres- 
ence of K-K' scatte ring; is determined by the range of 
scattering potential [20|. Short- (long-) range scatter- 
ers do (not) see the difference between A and B sites, 
and also, do (not) involve K-K' scattering. The weak 
localization behavior of graphene is indeed susceptible 
of presence or absence of K-K' scattering p3l . [2g]. In 
the absence of K-K' scattering, it is now established [23] 
that at one-loop order the graphene (a single Dirac cone) 
shows a weak anti- localization behavior, indicating that 
the system is symplectic (28| . 

An interesting question is to what extent this anti- 
localization tendency continues against generalizations. 
Recent numerical analyses suggest that this anti- 
localization tendency actually continues to the strong- 
coupling regime [24, 25], due probably to some cancella- 
tion of higher order terms in the expansion of I3{g). Ref. 
[2^ argues that such cancellation of higher order terms 
is a consequence of a non-trivial spectral flow, associated 
with Z2 Kramers symmetry, from which they coined the 
word, "Z2 (topological) metal". It is also pointed out 
that unconventional behaviors of Z2 metal can be also 
casted in terms of a topological term in the effective a- 
model description [29]. I n this paper, we achieve another 
generalization of Ref. [23] by reinterpreting the weak anti- 
localization property of graphene as the massless limit of 
a more general system, i.e., that of the Kane-Mele quan- 
tized spin Hall (QSH) insulator. 

Let us go back to the explicit form of Eq. ^ , and fo- 
cus on its properties under time-reversal operation. The 
time reversal symmetry (TRS) plays an essential role in 
the discussion of symmetry classes in weak localization 
theory. In the absence of K-K' scattering, the two Dirac 
cones are decoupled dynamically, whereas, TRS opera- 
tion transforms K to K' ^ and vice versa. In such cases, it 
is convenient to introduce the idea of "pseudo time rever- 
sal symmetry" (PTRS) [s^, in which one pretends that 
a represents a real spin so that pseudo-Kramers' pairs 
are formed in a single Dirac cone. Then PTRS trans- 
forms p to —p within a valley. Since r is invariant under 
PTRS by definition. Hi is also invariant. Unlike the real 
spin, however, the pseudo-spin should have single- valued 
eigenstates. This apparent difficulty is resolved in terms 
of the Berry phase [20.] as discussed shortly toward the 
end of this section. 



1. Graphene and its localization properties 

Graphene, an isolated single layer of graphite, has a 
band structure, with two massless points (often referred 
to K and K') in the first Brillouin zone, and in the vicin- 
ity of these points the low-energy effective Hamiltonian 
reduces to a Dirac- Weyl form [1]. This part of the Hamil- 
tonian, i.e.. Hi in Eq. ([1]), comes from the standard 



2. Topological mass and the quantized spin Hall effect 

In the second row of TABLE I shows, the two Dirac 
cones at K and K' valleys of graphene acquire a gap 
(the total Hamiltonian becomes H = Hi -\- H2) either in 
the presence of (a) imaginary hopping between second- 
nearest neighbors due to spin-orbit interaction [2] , or (b) 
AB sublattice symmetry breaking staggered chemical po- 
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TABLE I: Three strory structure of the Kand-Mele Z2 topological insulator, In addition to the real spin s, there appear two 
types of pseudo-spins: a representing A-B sub- lattices, and r specifying the valley: K oi K' . (*) If Fermi level is in the gap 





long-range scatterers 
(single valley) 


short-range scatterers 
{K-K' coupled by inter-valley scattering) 


(i) ideal graphene 
(massless) 


Hi = hVpiPxCTx +PyCry) 

single Dirac cone 


Hi = hvpiPxCrxTz +Pyay) 
opposite chiralities at K and K' 


(ii) mass terms: (a) topological, (b) ionic 
(a) QSH, (b) ordinary insulator (*) 


H2 = maz 


(a) H2 - -/^(TzTzSz = i/A 
(b) H2 = M(Tz = Hm 


(iii) Rashba spin-orbit interaction 
Z2 topological insulator 


Hr — —\R{(JxTzSy — GySx) 

mixes real spin | and [ (spin rotational- symmetry broken) 



tential. In case (a), the so-called topological mass term 
H/\ is generated, whereas (b) leads to a standard ionic 
mass term Hm- 

First note that in case (a), the topological mass term 
i^A = —^o-zTzSz is time reversal invariant. This stems 
from the fact that spin-orbit interaction preserves TRS. 
Second, note also that this is true only when we count 
both the (real) spin up and down sectors. Namely, if we 
pick up, say, only up spin part, then the system breaks 
TRS, showing, e.g., a finite axy (= ±e^//i, see below) 
even in the absence of magnetic field [33]. On the other 
hand, (b) a staggered chemical potential generates a mass 
term of the form: Hm = Maz, which has the same sign at 
K and K' points. The total Hamiltonian H = Hi -\- Hm 
describes an ordinary insulator such as a monolayer of 
boron nitride (BN). 

Consideration on such different types of energy gap 
(3^ naturally leads to the idea of "quantized" spin Hall 
insulator [2]. In contrast to the ordinary ionic mass case 
(b), where the so-called parity anomaly cancels between 
the K- and i^'-points and does not manifest itself [32], 
the topological mass term (a) has an opposite sign at 
K- and i^'-points [s^ as well as for up and down spins, 
resulting in the quantized spin Hall effect. There are 
actually two copies of quantum Hall states a la Ref. [s^ 
under zero magnetic field, one with spin t and cr^ = 
+e^//i, and the other with spin | and a^y = —e^/h [2]. 

The existence of an energy gap is also suggested in 
experiments. Its magnitude is under debate in photoe- 
mission experimen ts (341 , |35] In contrast to the theoretical 
prediction [23, 24, [23, weak localization experiments on 
graphene [3y, [33, [slf show also a unitary behavior. It 
is, therefore, natural to ask how the localization proper- 
ties would be influenced by the presence of mass term. 
Absence of WL may be attributed to ripples [sl, |^ . 

The type of mass term, given either by H/^ or Hm, is 
a relevant factor in our discussion on localization prop- 
erties (see later sections). If there is no K-K' scattering, 
and two Dirac cones are decoupled, however, the system 
cannot see the difference between the two types of mass 
term. The behavior of the system is thus strongly de- 
pendent on the presence or absence of K-K' scattering. 



We will see in later sections, when K-K' scattering is 
switched on, the topological mass (b) leaves the system 
unitary, whereas the ordinary mass drives the system to 
orthogonal symmetry class. Our analysis will be summa- 
rized in the language of PTRS in the final section (see 
TABLE II), and should be applicable to two dimiensional 
Z2 insulators in general. 



3. Rashba spin-orbit interaction and Z2 topological order 

Let us finally consider the third row (iii) of TABLE I, 
in which the total Hamiltonian becomes H = Hi -\- H2 -\- 
Hr. The Rashba term Hr turns out to be a relevant 
perturbation to the above symmetry properties. First, 
as for topological properties of the undoped phase, the 
quantized spin Hall effect is not robust, but replaced by 
a Z2 topological order. Q 

As for weak localization properties of the doped phase, 
the Rashba term Hr^ fixing the relative angle between 
Sz and the real space coordinates, changes the symme- 
try class. We will see in later sections that the Rashba 
spin-orbit interaction turns the system from unitary to 
orthogonal in the absence of K-K' scattering, whereas in 
the presence of K-K' scattering, the system turns from 
unitary to symplectic with weak anti-localization behav- 
ior. In the present "poor man's" analysis, we can identify 
the scattering channel that separates these different sym- 
metry classes. The two weak anti-localization phases, one 
in the graphene limit (single Dirac cone, unconventional) 
and the other associated with a Z2 topological insulator, 
evolve from each other via either orthogonal or unitary 
behavior, activated by the Rashba term together with 
K-K' scattering or the topological mass, respectively. 

Let us emphasize that the weak anti-localization be- 
havior of graphene occurs in the phase diagram of no 
K-K' scattering, wheras the weak anti-localization of Z2 
topological insulator occurs due to K-K' scattering. In 
the previous numerical analysis [l^, [l9| , this important 
fact has not been noticed because the disorder in the 
real-space model always involved the effective inter- valley 
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scattering. 



For example, we will use notations such as, 



\Ka) = \p+) 



COS ■ 



B. Construction of eigenstates 

To construct eigenstates explicitly, we first consider the 
simplest non-trivial case with vanishing Rashba interac- 
tion, and allow only long-range scatterers. Since A^^ = 0, 
real spin up and down sectors become decoupled. In 
terms of the Hamiltonian, H = Hi -\- H2 is diagonal in 
real spin 5-space as well as in r-spin space. We can, there- 
fore, consider separately = 1 (i^- valley) and = — 1 
(i^'-valley) . Taking = = 1 for simplicity, one may 
rewrite the Hamiltonian in the i^- valley, H = Hi -\- H/\ 
(with Tz = 1, Sz = 1), in the following simple form: 



H 



(2) 



by introducing a fictitious three-dimensional momentum 
P = {PxiPy, ~^)' As for H = Hi -\- Hm^ one has simply 
to replace it with p= {px^Py^ M). 

The Hamiltonian can then be diagonalized by choosing 
a proper quantization axis of the pseudo-spin, in analogy 
to the SU(2) spin case. One must take into account here 
that the a represents only a pseudo spin, as it is derived 
from the sublattice degree of freedom. Correspondingly, 
the momentum p, specifying the quantization axis for a 
is single-valued. Taking this single-valuedness also into 
account, one may denote the eigenvalues and the eigen- 
vectors of H = Hi -\- H2 as 



sin2 
-e*"^ cos ■ 



(3) 



where 6 and (j) are polar angles satifying 

-A , px 



cost/ = 



pI^pI^a^ 



I , cos (p - 



pI^pI 



(4) 



Here corresponds to the upper- (lower-) band eigen- 
vector. In the course of an adiabatic evolution of \p(t)±) 
around the origin of p, however, a Berry phase tt enters 
per winding. [41] This situation keeps consistency with 
the double- valued SU(2) eigenstates of a real spin [42]. 

In the following we will focus on the upper (conduc- 
tion) band with E = \p\ > A. One can also write 
down the eigenstates in the i^'-valley using the same 
parametrization. Because = — 1 in the i^'-valley, the 
conduction band has eigenfunctions of the form \p—) in 
Eq.([3|) with (j) replaced by —cj). We also introduce the 
notation, a, a' ^ etc. to specify the momentum p , in 
order to keep the consistency in notation with Ref . [23l| . 



\ — cos 



(5) 
(6) 



As long as the Rashba term is absent, one can safely fix 
s, say, to be j. Then we do not explicitly consider the 
real spin until Section HVl 



III. WEAK LOCALIZATION PROPERTY — 
UNITARY CASES 

Weak localization (WL) phenomena have been known 
since three decades [43]. Scaling to metal (weak anti- 
localization) was shown to be possible due to scattering 
by spin-orbit interaction [i^]. Absence of WL (unitary 
behavior) is, on the other hand, attributed to explicit 
breaking of time reversal symmetry (TRS). This paper 
shows that, in systems of graphene, a zoo of such different 
localization behaviors appears under the same Hamilto- 
nian, simply by activating or inactivating effective spin 
degrees of freedom (FIG. [1]). Graphene thus provides 
a contemporary aspect to the conventional WL theory 
framework. 

single Dirac 
(massless) 

AL 




FIG. 1: (Color online) Weak localization phase diagram of the 
doped Kane-Mele model in the presence of potential scatter- 
ers. Relation to other graphene-based models with different 
types of the mass, or ripples is taken into account. WL, AL 
and U refer to weak localization (orthogonal class) , weak anti- 
localization (symplectic class) and absence of WL (unitary 
class), respectively. 

We apply standard diagrammatic techniques for weak 
localization to the doped Kane-Mele model. In real sys- 
tems, such as a graphene sheet, doping can be easily 
done by simply applying a gate voltage. Suppose that 
our Fermi level is in the conduction band, so that the 
system is metallic in the clean limit. We then introduce 
weak disorder by taking into account scattering by impu- 
rities. To characterize our two-dimensional system of size 
L^, we focus on its longitudinal conductivity axx = g{L). 
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Weak localization refers iol/g correction to the so-called 
scaling function, 



d\ogg 
dlogL 



D 



2-^ 



(7) 



where /3{g) ^ {g ^ oo) in two spatial dimensions 
{D = 2). 1/^-correction to P{g) leads to logarithmic cor- 
rections to conductivity. Thus, the sign of 1/^-correction 
in Eq. ([7j), i.e., the sign of ci determines system's weak 
localization property: (i) Ci > 0: WL (orthogonal), (ii) 
ci = 0: absence of WL peak in resistance data (unitary), 
and (iii) ci < (symplectic). Symmetry classes (orthog- 
onal, unitary, and symplectic) are due to classification 
of corresponding random matrices [2^. Such logarith- 
mic corrections to conductivity can be calculated using 
diagram techniques based on Kubo formula [45]. Dia- 
grams contributing to a weak localization correction are 
particle-particle type ladders, or sometimes also called a 
"Cooperon". In contrast to the particle- hole type dia- 
grams, which always show a "diffusion type" sin- 
gularity, whether Cooperon diagrams are susceptible of 
such singularity is a more subtle issue related to 

time reversal symmetry of the system [43]. 

In this section, we first switch off Rashba spin-orbit 
interaction, and study whether different types of mass 
term lead to the absence of WL peak (unitary behavior). 
We consider only scalar potential scatterers in contrast to 
Ref. Q. Instead, we distinguish impurities of different 
potential range, and classify them into two categories, 
depending on whether they involve inter-valley scatter- 
ing or not. Long (short )-range scatterers involve (do not 
involve) K-K' scattering. 




limit 
AL (symplectic) 

^ E 

unitary 
\^NL (orthogonal)^ip 

^R//im ^ 

<^E = A-Xr ^ 

FIG. 2: (Color online) Weak localization properties for long- 
range scatterers (LRS) without K-K' scattering. The or- 
dinate shows crossover from weak anti-localization (WAL) 
to weak localization (WL) tendency {Xr = 0, as E is de- 
creased from the graphene limit toward the bottom of the 
band). Crossover from unitary to orthogonal symmetry class 



turns the system from weak localization to weak anti- 
localization. In the present case, its complex nature does 
not come from the scattering potential, but from the 
property of wave function. 

From the self-energy diagram, we define the scattering 
time tl' 



1 

TL 



27ru{E)nLul{\{Kp\l\Ka)\^) 



VL I cos"^ - + sin"^ - 



(9) 



A. Long-range scatterers 

Long-range scatterers involve only intra- valley scatter- 
ing. One can, therefore, safely focus on, say, the i^T- valley. 
Such scatterers do not distinguish between A and B sub- 
lattices, either, i.e., the impurity vertex is also unity in 
the AB sublattice space. The matrix element associated 
with long-range scatterers is, therefore, proportional to. 



{Kf3\l\Ka) = cos^ - + e^^^--^/^) sin^ 



(8) 



where we assumed elastic scattering {6^ = Op = 0). Note 
that here we pick up only the spin part of the matrix 
element, and \Ka) represents only the spinor part of the 
electron wave function. In order to find the full matrix 
element between different momentum eigenstates, Eq.(|8]) 
should be appended by by a spatial part that requires 
the momentum conservation. The phase factor in Eq. 
(j8|) is analogous to the Berry phase, which has already 
appeared in Ref. in the graphene limit {0 7r/2). 

The Berry phase, in the presence of a mass term, is not 
generally tt in contrast to the massless limit. Note also 
that Eq.(|8]) involves an imaginary number that indeed 



where n^, ul and i^{E) are, respectively, the impurity 
density, the strength of impurity scattering potential, and 
the density of states at the given energy E. They com- 
bine to give r]L = 2t:v {E)nLu\^ and (• • • ) represents the 
angular part of impurity average. 

The transport relaxation time Ttr involves the cos((/)q, — 
0/3) term as 



— = r]L{\{K(3\l\Ka)\^[l-cos{^a-^^)]) 

'Ttr 



VL 



COS — h sm - 
2 2 



(10) 



where we have used (cos((/)q, — (/)/3)) = 0, (cos^((/) 0,-^/3)) = 
1/2. Usually, the factor 1 — cos {(j) a — 4^ f3) in the first line of 
Eq. (p!Qj) is inactive for ^-function like (5- wave) scatterers. 
Here, the matrix (Dirac) nature of the Hamiltonian in- 
duces a cosine term that corresponds to vertex correction 
in diagrammatic language. 

The bare vertex function 7 has also (/)-dependence: 



4 ^ 
cos - 



sm 



-0/3) sin^^ 



= 7(0) 



.^(i)e^(0a 



-0/3) 



.^(2)g2i(0« 



-0/3) 



2 

(11) 
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In the second line, we classify terms according to their 
relative angular momentum, i.e., different (/)-dependence: 
where / = 0,1,2. This expansion helps to 
solve the Bethe-Salpeter equation: 



tanl9/2: 



■ af3 



(12) 



To find the solution of Eq. ([T2]), we expand also T into 
different angular momentum components: 



{<i>c-<t>,3) 



(13) 



and integrate over 0^, i.e., over the intermediate angle 
dependence in Bethe-Salpeter equation. Notice also that 
n :^ rL(l — TLDq^)^ with D being the diffusion constant: 
D = one finds, at the same order of precision. 
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(0 



(14) 



The crucial issue is the cancellation (or not) of the lead- 
ing order (~ 1) term in the coefficient of T^^\ In the 
graphene limit {6 = 7r/2), F^^^ becomes ~ l/g'^-singular 
driving the system to weak anti-localization. At the bot- 
tom of conduction band, on the other hand, F^^^ becomes 
more important, and toward the limit ^ ^ (though the 
model becomes ill-defined in this limit), it tends to show 
~ 1/g'^-singularity, leading the system to weak localiza- 
tion. Away from these limits, the system shows a unitary 
behavior, i.e., all the three Cooperons acquire a finite life- 
time. 

Here, let us recall that the weak localization refers to 
a logarithmic correction to the longitudinal conductivity 
axx, of the form. 



^cr^x ^ TA log 



(15) 



(L: size of the system, I: mean free path. A: constant 
of order /h). In front of the logarithmic term, — sign 
should be chosen for the weak localization case. Whereas, 
in the case of weak anti-localization, this overall sign in 
front of the logarithmic divergence is positive (the correc- 
tion tends to increase the conductivity). In the unitary 
case we mentioned above, the logarithmic divergence of 
Eq. (p!5|) as L ^ oo is cut off by the longest lifetime 
of a Cooperon. And in that sense, the behavior of the 
system is driven by the Cooperon of the longest lifetime. 
This situation is analogous to the case of spin-dependent 
scattering studied in Refs. [i^. In order to quantify 
such unitary behaviors, we define the lifetime r^^^ of a 
Cooperon, F^^^ such that 



1 



(16) 



where all r^^^'s are written in terms of a single parameter 



r(0) 
r(i) 
r(2) 



tan^ -, 



1 / 1 e_ 

2 \iM0m~^^''2 



tan4(l9/2)' 



(17) 



Since we have tan(^/2) < 1, there is no chance for F2 to 
dominate. The parameter 6 is transformed to energy E = 
A/ cos^. As E increases from E = which means that 
increases from 6> = 0, a crossover occurs at "universal" 
value. 



2 arctan —= 
V2 



1.23 



(18) 



When Oc < < 7r/2, F^^^ (weak anti- localization) be- 
comes dominant , leading to a positive logarithmic cor- 
rection to the longitudinal conductivity, of the form of 
Eq. (p!5|) with the + sign, but L replaced by t^^\ When 
6 < 60 F^^^ (weak localization) is dominant, and the cor- 
rection to cFxx is given by Eq. ([15]) with the — sign. In 
terms of E^ the crossover occurs dit E = Ec = 3A, and 
at this point the logarithmic correction changes its sign. 
This crossover behavior is illustrated in FIG. 1, where we 
consider A^^ = for the moment. 

Since long-range scatterers do not involve K-K' scat- 
tering, and two Dirac cones are decoupled, the system 
cannot see the difference between the two types of mass 
terms. Thus for both types of mass terms, we found 
predominantly a unitary behavior. Cooperons acquire 
a finite lifetime, which plays the role of cutting off the 
logarithmic correction, Eq. ([15]). Typically, there is 
no l/^f-correction to the /3(^)-function, in the L ^ 00 
limit. However, one can still see a crossover from anti- 
localization tendency to weak localization regime, as far 
as the longest lifetime of a Cooperon is much larger than 
the system size L. When Fermi level is far above the gap, 
the electron does not feel very much that there is a gap, 
i.e., he has a tendency to behave as if he were a mass- 
less Dirac fermion (weak anti-localization behavior in the 
graphene limit, E ^ 00). As the Fermi level approaches 
the bottom of the band, the electron starts to feel the 
gap, and when he is close, he even forgets about that he 
is actually a "relativistic fermion", and starts to behave 
(say, E < 3A) as if he were a non-relativistic electron, 
showing the weak localization behavior. 



B. Short-range scatterers with K-K' scattering 

Short-range scatterers couple different valleys, i.e., K 
and K' . One may, therefore, possibly see the difference 
between two different types of masses; ordinary and topo- 
logical. The scattering matrix elements involve also a 
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projection operator in the AB sublattice space, i.e., 



Matrix elements of such projection operators in the K- 
valley are, 



{KP\rA\Ka) = cos^ 



{Kf3\pB\Ka) = e 



(20) 



Matrix elements of such projection operators involving 
i^T'- valley depends, on the contrary, on the type of mass. 



7c 



It 



K Va K. 

Pa 
Pa' 



i PP 
I Pp' 



K' Va K' 
Vat- K' 



Pa 
Pa' 



i Pa 

{nAu\ 
_!_J^' 



Pa 
Pa' 



K' 



, K 
^ PP 



Pc 
Pa 



K' Vat+ K K' 



K' 
K' 



"trans" ^ 



FIG. 3: Particle-particle ladders. Bare and dressed Cooper- 
ons. Relevant diagrams in the KK' sector, "cis" and "trans" 
refers to specific configurations of the valleys: K and K' . 



1. Topological mass case 

Let us first consider the case of Kane-Mele quantized 
spin Hall insulator, for which the conduction band eigen- 
kets are given by Eqs. (|5|6p . Since we are concerned 
about K-K' scattering, let us first consider the inter- 
valley scattering matrix elements: 

{K'l3\VAr-\Ka) = e"^^/^ sin ^ cos ^, 

{K(3'\VAr+\K'a') = e^^«' sin ^ cos ^, 

{K'p\VBr-\Ka) = -e^^^ sin ^ cos ^, 

{Kp'\VBr+\K'a') = -e"^^"' sin ^ cos ^, (21) 

where r± = {tx ^iry)/2 are "spin- flip" operators associ- 
ated with the valley-spin {K-K'). Their contribution to 
scattering time is, e.g., 

27ru{E)nAu\\{K'(3\VAr-\Ka)\^ 
+ 27ru{E)nBul\{K'p\VBr-\Ka)f 

= 275sin2^cos2^, (22) 

where we have defined 

rjs = 2n.iE)^^^<±^. (23) 

nA,B and ua,b are, respectively, the impurity density 
and the typical strength of scattering potential at the 
A (B) sites. In order to obtain the full expression for 
scattering time, one has to consider also the contributions 
from intra-valley scattering, such as, 

27n. {nAul\{K/3\rA\Ka)\^ +nBul\{K/3\rB\Ka)f) 

( 6 6\ 

= 27rz/ I nAu\ cos^ — h tibu^ sin^ - ] (24) 



Here, we assumed that the strength of inter- valley scat- 
tering is the same as intra-valley scattering, but this sim- 
plification is irrelevant to our conclusions. The scattering 
time in the iC- valley reads, 

— = 27runAu\ {\{KP\VA\Ka)\^ + \{K' p\VAr-\Ka)\^) 

TK 

^2^vnBul {\{KP\VB\Ka)\^ + \{K' P\VBT-\Ka)\^){?h) 

and one finds a similar expression for tk' • One can verify, 
using Eqs. (|22|24p , that as far as nAu\ ~ ^bu% = ^su'^, 
one finds, 

I/tk I/tk' =r]s = 1/^5- (26) 

For short-range scatter ers we consider here, the trans- 
port relaxation time is identical to r^, since the projec- 
tion in the AB sublattice space leaves no cross term, i.e., 
0-dependent term; typically, ~ cos{<pa — 4^(3)1 in the ex- 
pression for 1/tk and I/tk'- 

As for particle-particle ladders, the momentum con- 
servation naturally leads us to classify them into KK^ 
i^i^'-mixed, and K' K' sectors. They correspond in the 
notation of Ref.[23], respectively, to J = 2, and —2 sec- 
tors. J = ja -\-ja' = j(3 -^j(3' is couscrvcd, where Jq, = ±1 
if a occurs in the K- {K'-) valley, since K — K' \^ only 
half of a reciprocal lattice vector. 

In the graphene limit, the i^i^'-mixed sector is 
most divergent. In the J = sector, two types of 
Cooperon diagrams are possible (see FIG. 2). Both of 
them have two K electron and two K' electron lines, but 
they appear either in the "cis" or "trans" arrangement 
(in the terminology of organic chemistry). Naturally, 7c 
(7t) refers to cis (trans), and the same rule applies to 
Tc^t' The key issue here is that as a result of projection 
Va,b and K-K' scattering, 7^ acquires an additional mi- 
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nus sign: 

7c = 2^vnAu\{Kp\VA\Ka){Kp'\VA\Ka') 
+ 27ivnBul {K'l3\VB\K'a) {K' I3'\Vb\K' a') 



-ft = 2nvnAu\{K'(3\VAT-\Ka){KI3'\VAT+\K'a') 
+ 27:vnBu\ {K' P\VbT- \Ka) {K[3'\Vbt+\K' a') 



— —e 



sin^(9 



(27) 



canceling with the Berry phase. Note also that as for 
0- dependence both 7c and have only the / = 1 com- 
ponent. Eq. ([27j) is a simple consequence of the matrix 
elements in Eq. ([2T|) . Recall also that (j)a — (t>a' = ^r, 
since + k^' = kp -\- kp^ = q c:^ 0. 

In order to calculate the correction to conductivity, 
we set /3 = a\ meaning k^ -\- kf3 = g :^ 0, therefore, 
— <I^I3 = TT. The Bethe-Salpeter equation becomes two 
coupled equations: 



(28) 

Indeed, both Tc and Ft contribute to the l/q^ singularity. 
After diagonahzation, one finds, 



Tc + rt = 0, 

(7W- 7^^)05! (Fe-r,) = 7c -7t, (29) 



where Us — rs{l — rsDq^). The cancellation of the 
leading order (~ 1) term is incomplete, giving a finite 
life time tkk' to the Cooperon: 




r, = 



D*(0«-0/3) 



2Tl{Dq^ + l/TKK')' 

The lifetime tkk' behaves, as a function of like 

TKK' -E2 - A2 



(30) 



(31) 



Clearly, this KK' (J = 0) Cooperon mode shows 
singularity only at the E ^ 00 limit [23]. Decreasing 
energy toward the bottom of the band, another KK (J = 
2)-mode becomes important. The lifetime of the latter 
Cooperon (KK-mode) behaves, as a function of E^ like 



Tkk 



E^ 



£;2 + A2 • 



(32) 



Thus, except at £^ = A and ^ ^ 00, all the Cooperons 
having a finite life time, the system shows a unitary be- 
havior with no 1/^-correction in the L ^ 00 limit (see 
FIG. 3). 




massless limit 
1 /tk k' WL (orthogonal) 



► E 



AL (symplectic) 



FIG. 4: (Color online) Topological mass case. Weak local- 
ization properties in the presence of short-range scatterers: 
inter- valley scattering is allowed {K-K' coupled). 



One may find this contradictory to the fact that the 
topological mass term — Aa^r^s^ does preserve TRS. To 
clarify this point, first note that the entire Hamiltonian, 
Eq.([T]) is time reversal invariant, and that this comes 
from the fact that spin-orbit interaction preserves TRS. 
Second, note also that this is true only when we count 
both the (real) spin up and down sectors^ and that if we 
pick up, say, only up spin part, then the system, show- 
ing, e.g., a finite Cxy (= /h) even in the absence of 
magnetic field [33], clearly breaks TRS. In the absence 
Rashba spin-orbit interaction, spin up and down sectors 
are actually decoupled, indicating that the system be- 
longs to the unitary class. 



2. Ordinary mass case 

In the case of ordinary mass, one has to replace \K'a) 
in Eq.(|6|) with 

\K'a) = . (33) 

In this basis, inter-vahey matrix elements become 

{K'f3\pAT-\Ka) = e-^'^-'cos^^, 
{Kp'\VAT+\K'a') = e'*"'cos2^, 



{K'l3\VBT-\Ka) = -e*Wsin' 
{Kp'\VBT+\K'a') = -e-'*='sin2^. (34) 

Notice that their contribution to the scattering time is 
identical to Eq. ([24|) . One finds, therefore, substituting 
Eqs. ([20|) and (|3l|) into Eq. ([25|), that the scattering 
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time becomes this time, 



— = An v {E) ( uaUj^ cos --{-ubu^ sin - 
tk \ ^ 2 



2r]s cos - +sin - = — . 



(35) 



One can also verify that 1/tk' = On the other 

hand, the expressions for 7c, t become also. 



7c,t = ±e*(^«-^^) 



2 4^1 2 • 4 ^ 

uau ^cos -^UBU^sm - 



X nu% 



4^1 • 4 ^ 

COS — h sm - 
2 2 



(36) 



Notice that Eqs. (|28|29p are always valid, whereas here 
IIc^t's{1- r'sDq^). Eqs. ([35l ^ suggest that in con- 
trast to the topological mass case, the scattering time 
(the self energy) and the bare vertex function (times 11) 
always cancel identically (giving unity) at the leading 
order of Eq. ([29|) . As a consequence, the Cooperon dia- 
gram shows singularity, which occurs always at the 
/ = 1-channel: 



Ft 



(37) 



Tt is indeed positive and singular, independent of 9^ in- 
dicating weak localization whenever the Fermi level is 
above the gap. 

We have thus seen a clear distinction between two 
types of mass terms in their weak localization proper- 
ties, as K-K' scattering is switched on. In the ordinary 
mass case, the system shows the orthogonal behavior, 
in sharp contrast to the unitary behavior of topological 
mass case. Such different weak localization properties 
due to different types of mass terms can be understood 
as follows. Recall that the topological mass term induces 
a finite (Jxy^ if one picks up only one of the two real spin 
components. That was a clear signature of broken TRS, 
leading to unitary behavior. In the ordinary mass, on the 
other hand, contributions to Gxy from each valley cancel, 
i.e., (j^y + (j^y = 0. Therefore, there remains no trace of 
broken TRS any more, once two valleys are coupled by 
K-K' scattering. The ordinary mass term does preserve 
TRS, and as a result one finds always a diffusion- type 
singularity, leading to a weak localization behavior. 



IV. RASHBA SPIN-ORBIT INTERACTION 

Rashba spin-orbit interaction Hr is also an impor- 
tant factor for characterizing the physical properties of 
the doped Kane-Mele model. A finite Rashba term ap- 
pears only when the system loses its inversion symme- 
try along the z-axis (perpendicular to the 2D "graphene 
layer"). Physically, such breaking of inversion symmetry 
can be introduced by an asymmetric potential, e.g., when 



a graphene sheet is placed on a substrate. If Rashba 
spin-orbit interaction is stronger than a critical value 
{\r > 2A), the system actually has no topologically non- 
trivial phase \3\ . Here we suppose that Rashba spin-orbit 
interaction is not too strong, and the undoped system is 
still in the topologically non- trivial phase. However, we 
show, in this section, that in regard to weak localization 
properties of our doped system, Rashba spin-orbit inter- 
action is still a relevant perturbation, and changes the 
symmetry class of our system, as soon as it is turned on 
(in real samples of finite size of order L^, there will be a 
crossover at the strength of Rashba spin-orbit interaction 
of order A^^ ~ 1/^)- This is because the Rashba spin- 
orbit interaction mixes real spin up and down. Before 
switching on A^^ 7^ 0, we did have real spin up and down, 
but they were just there, and inactive. 

In the presence of Hji, we can still work on a 4 x 4- (in- 
stead of 8 X 8-) matrix space associated with a and 5 since 
Tz is diagonal in our basis. Focusing on one of the two 
valleys, say, K, one notices that the Rashba spin-orbit in- 
teraction couples A| and B| only (here, the up and down 
arrows refer to the real spin). The total Hamiltonian in 
the iC- valley reads. 



Hk 





-A 


Px - iPy 








Px 


^iPy 


A 














A 


Px - iPy 


\ 








Px + iPy 


-A 



(38) 



where the inner 2x2 structure refers to AB-spin, whereas 
the outer 2x2 block structure is asoociated with the real 
spin. Here, for the sake of simplicity, we consider only the 
case of Kane-Mele type topological mass term: AazTzSz. 
After diagonalization, the resulting four energy bands are 
classified into two conduction, and two valence bands, 
which we will call, respectively, and The two 

valence bands E(i± are degenerate on their top: the top 
position is always at £^ = — A, unaffected by the Rashba 
term, whereas the conduction bands Eu± are split by 
2A,^: 

Eu± = \p±\±Xr/2 



= ^P^+pg + (A±A^/2)2±A^/2, 
Ed± = -\p±\±Xr/2 



-Jpl^pl^{A±Xn/2y±Xn/2. (39) 



The Corresponding eigenvectors \u^) and can be 

conveniently parametrized by introducing a fictitious 3D 
momentum: 



P± 




(40) 
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and associated polar angles: 



cos 9± 



A ± Afl/2 



Px 



pI+pI 



(41) 



Note that (j) is actually common to all cases. Thanks to 
the parameters introduced in Eq. (jHj), the eigenvectors 
corresponding to the eigenvalue Eu^d±i given in Eqs. ([39|) 
allow for the following compact representation: 



\Ku±) 



\Kd±) 



1 

71 



1 

71 



"sin% \ 



COS ■ 



cos -| 
y ^ie^^ sin - 



'*cos% \ 
■sm^ 



±1 sm 



\ Tie"*' cos % / 



(42) 



Note that a crucial difference here compared with the 
unitary limit (A^^ = 0) is that because of the Rashba 
coupling, we have a stronger constraint on the choice of 
our basis, and as a results there is no Berry phase in 
the matrix elements (see below). In the following, we 
consider again a doped system (our Fermi level is in the 
conduction band), i.e., E > A — Xr. 



A. Long-range scatterers: weak localization 

The long-range scattering impurities do not couple K 
and valleys, and in a given valley do not distinguish 
A and B sites, i.e., its scattering potential is propotional 
to unity (no projection). In the following we focus on the 
K-valley. 



1. A-Xr<E<A + Xr case 

Let us first consider the case in which the Fermi level 
is close to the bottom of the conduction band, A — A^^ < 
^ < A + Ai^. In this case only )-branch contributes 
to our diagrams. With this remark, we will omit, in the 
following, the indices u and — , for specifying a ket such as 
\Ku — a)^ and denote it simply as \Ka). As we have seen 
in Eqs. (|42|) . because of the Rashba coupling, we have 
now a stronger constraint on the choice of our basis, and 
as a results there is no Berry phase in the matrix element: 



given in Eq. (|4Q|). and specified by a. Self-energy dia- 
grams give the bare scattering time (without vertex cor- 
rection) as 



1 . 4 

-sm 



cos 



(44) 



Because of the cos{<pa — 4'i3) term in Eq. (|43|) the bare (and 
also dressed) vertex function has several angular momen- 
tum components. Naturally, we expand them such that 



la(3 



■ a(3 



(0 

(3^ •> 



(0 Ji( 



(45) 



The bare 7/'s are given as. 



v(0) - 



v(l) 



1 4 
-COS 



sm 



.(-1) = 'JL 



— COS i 

4 



7 



(2) 



7 



(-2) VL . 4^- 

^ > — — sm — - 
4 



(46) 



To solve the Bethe-Salpeter equations, first recall that 
there is an angular integration over the intermediate an- 
gle 0^, which forbids coupling between F^^^'s with differ- 
ent angular momentum I. The Bethe-Salpeter equations 
take the form of Eq.([Tlj). Then, comparing Eqs. (jHj) and 
(j46j) , one can verify that Fq shows a diffusion- type 1 / c^- 
singularity, with a positive amplitude, leading to weak 
localization. 

Recall that in the absence of Rashba spin-orbit interac- 
tion, the anti-localization tendency toward the graphene 
limit was given by the / = 1 term. We emphasize 
that here the relevant (singular) contribution is from the 
/ = 0-channel, and this is quite contrasting to the former 
case. We mentioned earlier that as Rashba spin-orbit in- 
teraction Ai^ 7^ 0, with off-diagonal matrix elements in 
the real spin space, is turned on, the Berry phase, as- 
sociated with the Xjq^ singularity of the / = 1-channel, 
disappears. What we have discovered above is consistent 
with this observation. It is also quite natural (as far as 
K'- valley is switched off) that the Xjc^ singularity ap- 
pearing in the / =(even) channel leads to weak localiza- 
tion, whereas in the / =(odd) case, the same singularity 
leads to weak anti-localization. From a perspective point 
of view, it is also useful to remark that the number of 
effective spin degrees of freedom activated in the system 
is increased from A/'g = 1 (A^^ = 0) to A/'g = 2 A^^ 7^ 0, 
in switching on the Rashba term. This aspect is summa- 
rized in TABLE II as explained later. The parity of 
is a decisive factor for determining the symmetry class of 
system. 



(g) \\Kol) = cos(0c. - ^/3) sin^ — + cos^ — , (43) 

which is indeed real. Note that \Ka) here denotes \Ku—)^ 
given in Eqs. (j32j), for a fictitious 3D momentum pi. 



2. A + Xr < E case 

When the Fermi level is above the bottom of upper 
branch of the conduction band, i.e., when A -\- Xr < E, 
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Pc 
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^PP' 
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FIG. 5: Particle-particle ladders involving inter-branch pro- 
cesses. In the presence of finite Rashba term Xr ^ (a) 7 — y 
does contributes to singularity, whereas such diagrams 
as (b) and (c) are irrelevant to the singularity, since + Pa.' 
cannot be smaller than the order of A^. 



both \Ku—) and \Ku^) channels contribute to the weak 
localization properties. In order to shorten the equations, 
we omit in this section even the index K in the brackets, 
and keep only ± for specifying branches, e.g., | — a) = 
\Ku — a). In this new notation. 



{-p\l\-a) = cos(0c. - (^/3)sin^ — + cos^ — , 



= -isin((/)c-' 
= cos((/)a - 0/3~ 



(7 1 U — 

) sm — sm — 
2 2 



sin2^ + cos2^.(47) 
2 2 ^ ^ 



Using these matrix elements, one can calculate the scat- 
tering time r± for the \Ku^) branch: 



1 



1 . 4 

- sm 
2 



1± 
2 



cos 



1± 

2 



- sm 
2 



■ sm 



(48) 

where the last term corresponds to the contribution from 
inter-branch matrix element, — a) p. 

The inter-branch matrix elements in Eqs. (|47|) also 
appear in the particle-particle ladders (see FIG. 4). In 
principle, the four electron states a, can be- 

long to either of the two channels, \Ku—) and \Ku-\-). 
There is, however, an important simplification at this 
level. Since we are interested only in the singular 
part of Cooperon diagrams, we need Pa + Pa' = g — 0, 
and similarly, pf^ -\- pp' = g — 0. This means that a and 
a' must belong to the same branch, which is similarly 
the case for P and /3^ We are thus led to consider such 

diagrams as, 7--, 7++, 7 and 7+-. 7-- has already 

appeared in the regime. A— Ai^ < E < A-\-Xr (its explicit 
form is also shown there, in the form of expansion with 
respect to 7 ^ is defined in FIG. 4 (a). 7^ is sim- 
ilar to 7 ^, only with \Ku—) and \Ku-\-) interchanged. 

Other 7's such as FIG. 4 (b) and (c) are irrelevant to 



singularity. Explicit form of 7's are given as. 



Vl 

2 COS' 



cos 

2/ J. 



I sm 



3) sin^ — cos^ h cos^ — 

^22 2 



7++ 



7-+ 
7+- 



VL 

2 COS' 



cos 

2/ J. 



4 ^+ 



I sm 



. 2 ^+ 2 ^+ 4 ^+ 

I sm — cos — + cos — 



-r]L sin^(0c. - (j)f3) sin sin — 
7-+- 



(49) 





The dressed Cooperon diagrams, r__ and F^ , satisfy 

a coupled Bethe-Salpeter equation, which takes the fol- 
lowing form: 



7__n_ 7-+n+ 
7+-n- 7++n+ 



where Il± :^ r±(l — D±T±q^). Recall that in the last 
term of Eq. ([5Q|) , combinations of the type, 7nr, appear, 
which implicitly contain averages over the 0-angle. We, 
therefore, expand bare 7's and dressed F's, into different 
angular momentum contributions, as Eqs. (|35|), and, to 
identify the singular contribution, pick up only the / = 
component. One way to convince oneself that the dressed 
Cooperons show indeed singularity at the / = 

channel, is to prove the following identity. 



det 



1 (0) (0) \ 

(0) (0) - ^• 



(51) 



In order to verify, first notice 



(0) 

7-- = 



1 4^- 
- COS — 

2 2 



sm 



(0) 

7-+ 



r]L . 2 ^- • 2 ^+ 

sm — sm — 

2 2 2 



(0) 

= 7;- 



(0) 

7;; = r]L 



-L 4 (7 + 

- COS — 

2 2 



sm 



4 ^+ 



(52) 



and use such relations as, 
(l/r_ - 



(0) 

■ 7__r_ 



(0) \ 



1 (0) 



r]L . 2 ^+ • 2^- 
T- — sm — sm — 
2 2 2 

(0) 

-^+71+- 



(0) 



(53) 



Based on these observations, we claim that Rashba 
spin-orbit interaction recovers the singularity of 

Cooperons, driving the system back to weak localization 
(orthogonal symmetry class) , whenever the Fermi level is 
above the gap. These features are illustrated in FIG. 1. 
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Let us finally consider what happens if one adiabati- 
cally switches Oj^ the Rashba term. The simplification we 
have made for justifying Eq. ([5Q|) is no longer valid. In 
the limit of vanishing A^^ we cannot simply neglect such 
diagrams as FIG. 4 (b), (c) and similar diagrams. In 
principle, they could contribute equally to the Xjc^ sin- 
gularity if the singularity ever appears. Relations such as 
Voi = can become satisfied in these diagrams. 

At the same time, inclusion of all such diagrams makes 
the size of coupled Bethe-Salpeter equation much bigger. 

One can no longer decouple r__ and as Eq. ([5Q|) . 

The Cooperons acquire more channels to couple to, and 
as a result, the cancellation property between particle- 
particle ladders and the self-energy, such as Eq. ([5T]) . is 
lost. This loss of cancellation property due to activation 
of those channels which are depicted in FIG. 4 (b) and (c) 
explains why the system becomes unitary in the A^^ = 
limit. 

On the other hand, it is also possible to choose a basis 
if Ai^ = 0, as we did in Sec. Ill, in such a way that 
two upper band branches are decoupled. Recall that in 
the absence of Rashba term, the real spin part of the 
Hamiltonian is diagonalized in the real spin basis (by 
diagonalizaing s^)- 



B. Short-range scatterers: weak anti-localization 

In the presence of Rashba spin-orbit interaction, we 
found weak localization behavior for long-range scatter- 
ers, on contrary to the graphene limit (compare the first 
(i) and third (iii) row of Table 1, left column). Note also 
that K-K' scattering drives the system, in the graphene 
limit, from symplectic to orthogonal symmetry class. We 
thus finally consider short-range scatterers in the pres- 
ence of Rashba interaction. Since short-range scatter- 
ers involve inter- valley scattering, we need to consider 
also the eigenstates at the i^'-valley. Diagonalizing the 
Hamiltonian at the iC'-valley, one veryfies that the en- 
ergy spectrum is identical to Eq. (|39|) : there are two va- 
lence bands degenerate on their top at £^ = — A, whereas 
two conduction bands split by 2Ai^. Using the same pa- 
rameterization as before, i.e., Eqs. (|4Q|) and (|4T]) . the 
corresponding eigenvectors read. 



\K'u±) = — 



1 



COS ■ 



2 

y ±i cos ^ J 



=Fie"^ sm • 



\K'd±) = 



1 

V2 



sm 



2 

-id ti- 

e ^ cos — 
iie*"^ cos 
V i^sin^ J 



(54) 



Short range scatterers do distinguish AB pseudo-spin, 
whereas as far as they are non-magnetic impurities, they 



are unconcerned about the real spin. The impurity po- 
tential operator is proportional to Va,b^ which should be 
understood here as. 




(55) 



Let us focus on the regime: A — A^^ < E < A + A^^, in 
which only \Ku—) and \K'u—) modes are available. Us- 
ing a, /3, • • • for specifying momenta, the relevant matrix 
elements read, 



{K'p\VAr-\Ka) 
{K(5'\VAT+\K'a') 



■ sm ( 



(56) 



As for the latter matrix element, we noticed (j)^' — ^a+^r, 
(/)^/ (/)^ + TT in the second expression. This clarifies the 
nature of additional minus sign analogous to Eq. ([3Q|) . 
As for particle-particle ladders, there are, as before, into 
KK^ KK'-mixed, and K' K' sectors. We focus again 
on the J = sector, in which two types of bare ver- 
tex functions, 7±, are possible, and correspondingly, two 
types of Cooperons: T± (recall FIG. 2). T± obey a cou- 
pled Bethe-Salpeter equation equation, ([5Q|) . Again, as 
a result of projection Va^T^b and K-K' scattering, 7^ 
acquires an additional minus sign: 



Vs 2 
— cos 
4 



■ sm 



(57) 



However, in the present case, the Berry phase had disap- 
peared, implying that this time the new minus sign leads 
the system from orthogonal to symplectic. To verify this, 
first notice that intra- valley matrix elements are given as. 



{K(3\VA\Ka) = 
{K'(3\VA\K'a) = 



-0/3) sin^ , 



cos 



2 

-0/3) 



COS 



sm 



These combine to give. 



7c = -r 



4 



— - +COs((/)c - (/)/3)- 



(58) 



(59) 



which is indeed positive. The scattering time is also cal- 
culated to be, 



1/tk 1/tk' Vs/^' 



(60) 



Note that all 6>_ -dependence cancelled, giving unity. 

Let us now sum up the particle-particle ladders, and 
then solve Bethe-Salpeter equation. After diagonaliza- 
tion, one finds, 



[1 - (7c - 7t)n] (Fe 



fTt ~ (regular), 
Tt) = 7c - 7* > 0. (61) 
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In order to identify the singular contribution, one should 
expand 7's and F's into different angular momentum 
components, and pick up the / = term. Other con- 
tributions are indeed regular. Comparing Eqs. (|6Q|) . ([59|) 
and ([57|) . one can indeed verify that Tc — Tt shows 
singularity, whereas only Ft contributes to the conduc- 
tivity, giving weak anti-localization correction. 

Based on this observation, together with our analysis 
in the long-range case for the regime: E > A -\- Xr^ one 
can naturally conjecture that this symplectic tendency 
continues to the higher energy regime: E > A -h A^^. 
Then, the phase diagram in the presence of K-K' scatter- 
ing (FIG. 3), becomes predominantly symplectic (when 
Ai? 7^ 0). Comparing two phase diagrams, in the absence 
(FIG. 1) and presence (FIG. 3) of K-K' scattering, one 
can see that the two weak anti- localization behaviors, one 
in the graphene limit, and the other of Z2-topological in- 
sulator phase, have a quite different origin. The former 
occurs in a single Dirac cone, whereas, the latter occurs 
due to K-K' scattering. The former is related to Berry 
phase a la Ref.[42^, whereas in the presence of Rashba 
term mixing real spins, the matrix elements become real 
(the Berry phase disappears). 

Let us finally estimate the strength of gate electric 
field required for observing the crossover to weak anti- 
localization. For the crossover to be experimentally ac- 
cessible, Rashba SOI needs to be the order of ~ 1 K. 
This corresponds to the electric field of order ~ IV/nm 
p, a value attainable in double-gated graphene devices 
[43]. The crossover to weak anti-localization will be ob- 
served for a sample with insignificant ripples. A similar 
crossover due to Rashba SOI has been observed in an- 
other context in InGaAs/InAlAs quantum well [48,]. 



V. SUMMARY AND DISCUSSIONS 

We have studied localization properties of the doped 
Kane-Mele model. We assumed that the disorder is weak 
enough to apply the standard weak localization theory, 
i.e., we have focused on the leading order 1/g correc- 
tion in the 1/g expansion of [3{g). We have consid- 
ered a phase diagram in the {E^ Ai^)-plane together with 
inverse lifetime of a Cooperon, 1/r, in order to indi- 
cate the crossover behavior in the unitary case. When 
r ^ 00, the corresponding Cooperon mode shows 
singularity and gives correction to the /3-function. Oth- 
erwise there is no correction to the /3-function at the 1/g 
level and in the L ^ 00 limit. In the presence of topolog- 
ical mass term, the system shows predominantly a uni- 
tary behavior, as expected from symmetry consideration. 
However, as we have seen in Sec. Ill, the role of mass 
term to system's localization property is a more subtle 
issue, depending on the type of impurities, and the pres- 
ence or absence of Rashba spin-orbit interaction. 

A broader perspective is obtained on the phase dia- 
gram of localization properties in terms of the number 
Ng of activated (pseudo) spin degrees of freedom. TA- 



BLE [ni summarizes the weak localization (WL) proper- 
ties of doped Kane-Mele model. It shows that contrasting 
behaviors are closely related to the parity of Ns. To be 
explicit, if Ns is even, the system shows WL, whereas if 
Ns is odd, the system undergoes weak anti-localization 
(AL). The unitary behavior, on the other hand, is inde- 
pendent of the parity of Ns^ and emerges whenever the 
effective time-reversal symmetry, i.e., TRS or PTRS, is 
broken. Thus, in order to identify the weak localization 
symmetry class, it is most important to count active spin 
degrees of freedom, and whether the effective time rever- 
sal symmetry is preserved. Note that the pseudo-spin 
a is odd under PTRS, but even under the genuin TRS. 
Therefore the ordinary mass term in the third column 
does not break the TRS, and the case (a) shows the or- 
thogonal behavior. 

In order to distinguish active and inactive spins, we 
introduce (pseudo) TRS operations T5], defined in the 
subspace S, such that 

Ts(ffi + Hr)T^^ = Hi+ Hr, (62) 

where S = {a}, {o^, r}, {(J, s}, {a, r, 5}. Their explicit 
representations are, 

T(j = icFyC^ Tfjj- = TxC ^ 

Tas = {-i(7y){-iSy)C, Terrs = rx{-iSy)C, (63) 

where C is complex conjugation. T^rrs represents the 
genuine TRS operation. Effective TRS of the system is, 
therefore, determined by the transformation property of 
the mass term (see TABLE lllip . When a mass term is 
odd against TRS, the system shows the unitary behavior. 
If some (pseudo or genuine) TRS exists in the system, its 
weak localization property is determined by the number 
Ns of the activated spin degrees of freedom. One can 
verify T| = 1 if Ns is even, whereas T| = —1 if Ns 
is odd. The former (latter) corresponds to the orthogo- 
nal (symplectic) class in the random matrix theory [28| , 
and leads to constructive (destructive) interference be- 
tween two scattering processes transformed from one to 
the other by Tj^. 

The above symmetry arguments allow for a slight gen- 
eralization of our localization phase diagram (see FIG. 
[1]). Recall that in Sec. IV A, we have chosen the mass 
to be H2 = —AazSz. Results of this analysis are con- 
sistent with the symmetry consideration. In FIG. [H we 
further conjecture based on the symmetry consideration 
that a similar analysis for an ionic mass term H2 = Ma^ 
leads to the absence of WL (unitary behavior). Thus, the 
four minus signs in TABLE IIIII correspond to four uni- 
tary phases in FIG.[T1 We also point out that in terms of 
symmetry, ripples plays the same role as the ionic mass. 

We mention that the results of our analyses summa- 
rized in TABLE [III and FIG. [1] are also relevant to in- 
trinsic single valley systems ^ such as the one realized 
in HgTe/CdTe quantum well Q. In the latter system, a 
single pair (Kramers pair) of Dirac cones appear, in con- 
trast to graphene (cf. graphene has two pairs of Dirac 
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TABLE II: Summary: weak localization properties of the Kane-Mele model. The weak localization (WL) or weak anti- 
localization (AL) results depending on the presence or absece of inter-valley scattering, and Rashba spin-orbit interaction. 
Notice the role of Ns, the number of activated (pseudo) spin species. The value of angular momentum / is also indicated for the 
singular 1/^^) Cooperon channel (see Eq. (|45|) ). In the presence of K-K' scattering, activation of the valley spin increases 
Ns by one, leading to change of weak localization properties. 





long-range scatterers 
(no K-K' scattering) 


short-range scatterers 

[K and K' points coupled) 


(i) graphene 


a: Ns^l (odd) 
1/^2 . / ^ 1 ^ AL 


(J, r: Ns =2 (even) 
1/^2 . / ^ 1 ^ WL 


(ii) finite mass system 

(a) ordinary 

(b) topological 


a: Ns^l (odd) 

unitary (WAL as ^ ^ 00) 


(J, r: Ns =2 (even) 

(a) : / = 1 ^ WL 

(b) unitary (WL as ^ ^ 00) 


(iii) doped Z2 insulator 


a, s\ Ns — 2 (even) 
1/^2 . / ^ ^ WL 


a, r, s\ Ns=3 (odd) 
1/^2 . / ^ ^ AL 



TABLE III: (Pseudo) time reversal operations Ts, relevant 
in the subspace spanned by activated spins. Transforma- 
tion property of a mass term O = mcr^, /\azTzSz under T^: 
T^OT-^ = ±0. The sign appears in the table. U refers to 
unitary class. 



activated spins 


a 




(J, s 


a,r, s 


relevant TRS operation 




T 

J- err 


T 

J- as 


T 


(Jz 




+ 




+ 


1\GzTzSz 






+ 


+ 



cones; they appear at K and K' points). Note that two 
Dirac cones (Kramers pair, with real spin tand |), de- 
scribed in the model introduced in Ref. [J], have the 
same sign for their mass term, implying (superficially) 
the ordinary (ionic) mass case in our language. However, 
they have also different chiralities, and by a simple linear 
transformation (exchange of rows and columns), they can 
be mapped, and indeed corresponds to the topological 
mass case, studied in Sec. IV A. Our diagnosis (see TA- 
BLE [III and also FIG. [T]) suggests that the single- valley 
HgTe/CdTe system shows a crossover from unitary to 
orthogonal (WL) symmetry class, on activating the real 
spin degrees of freedom by a (Rashba-like) off-diagonal 
interaction between real spin | and | sectors. 

What is then essential to the AL behavior character- 
istic to the doped Kane-Mele model? As can be guessed 
by the analysis relying on the number Ns of active spin 
species, the combination of Rashba interaction and the 
K-K' scattering makes A^^ = 3. In other words, the AL 
results even in the absence of the topological mass term, 
or even in the case of the ordinary mass term. In this 
sense, the AL is not a peculiar property of the doped Z2 
insulator. The hallmark of the doped Z2 insulator be- 
comes manifest in the absence of the Rashba interaction. 
Namely, the robustness of the unitary behavior against 



the range of the disorder potential is a fingerprint of topo- 
logical mass term. The unitary behavior is closely related 
to the quantized spin Hall effect in the undoped case. 

Finally, we comment on the role of the channel index 
/ which is a quantum number associated with a relative 
momentum of electrons before and after the scattering by 
an impurity. See Eqs. (pT]) and (|45]) for its definition. As 
we have already seen in the body of the paper, to identify 
the singular channel / greatly helps to determine the weak 
localization property. In the case of long-range scatterers, 
the two valleys are decoupled, and we can safely focus 
on, say, the i^-valley. Then, in Sec. Ill and IV, we have 
seen that the singularity appearing in the / =(even) 
channel leads to WL, whereas in the / =(odd) case, the 
same singularity leads to AL, This situation is reversed 
in the presence of K-K' scattering: even / leads to AL 
and odd / to WL according to TABLE II. 

The above selection rule on the singular channel / is 
superficially dependent on the choice of the spin part of 
the wave function. In Sec. II and III, we have chosen to 
use a single- valued two-component spinor for describing 
AB sub-lattice spin eigenstates, which were naturally ex- 
tended to the single- valued four-component spinor in the 
presence of Rashba spin-orbit interaction. If one choose 
a double- valued spinor instead of Eq.([3]), pretending that 
the AB-sublattice spin realizes a real spin, then the se- 
lection rule is shifted by one. This shift, however, is the 
same for all cases in TABLE II. Hence if we take the dif- 
ference of / relative to the graphene limit, for example, 
the WL or AL behavior does not depend on the choice 
of the spin eigenstates. 

As we have seen in Sees. Ill and IV, the value / of 
singular channel does not change in the presence of K- 
K' scattering. The value of Ns is, on the other hand, 
increased by one, in switching on the KK' valley spin 
r. Since r flip does not involve <pa variables, the value / 
of the singular channel remains the same. However, the 
nature of the channel does change by the K-K' scatter- 
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ing. Namely, switching on of the vahey spin appears as 
an additional minus sign in particle-particle ladders. 

In conclusion, we have constructed a poor man's phase 
diagram of the doped and disordered Kane-Mele model. 
In order to characterize the Z2 nature of the model, we 
considered a doped case in contrast to the more familiar 
topological insulator phase. We also characterized the Z2 
nature from its &i^/A; properties, instead of its eti^e proper- 
ties, to which the Z2 nature is often attributed. [2, 7] Our 
analysis is restricted to the weak (anti-) localization level. 



and can be used for constructing an effective a-model de- 
scription of the same model. The basic mechanism dis- 
cussed here will be also useful for the understanding of 
different types of topological insulators. 
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